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Abstract 



ergies, we study the production of dileptons with mass below about 1 GeV 



Using a relativistic transport model to describe the expansion of the fire- 
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\q , from these collisions. The initial hadron abundance and their momentum dis- 

o\ 

tributions in the fire-cylinder are determined by following the general features 
of the results from microscopic models based on the string dynamics and fur- 
ther requiring that the final proton and pion spectra and rapidity distributions 



are in agreement with available experimental data. For dilepton production, 



we include the Dalitz decay of 7r°, 77, rf, uj and a\ mesons, the direct decay 
of primary p°, uj and <p mesons, and the pion-pion annihilation that proceeds 
through the p° meson, the pion-rho annihilation that proceeds through the 
di meson, and the kaon-antikaon annihilation that proceeds through the (j> 
meson. We find that the modification of vector meson properties, especially 
the decrease of their mass due to the partial restoration of chiral symme- 
try, in hot and dense hadronic matter, provides a quantitative explanation 
of the recently observed enhancement of low-mass dileptons by the CERES 
collaboration in central S+Au collisions and by the HELIOS-3 collaboration 
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in central S+W collisions. 
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I. INTRODUCTION 



The study of the properties of vector mesons in hot dense matter formed in heavy-ion 
collisions is one of the most exciting problems in nuclear physics. Theoretically, lattice QCD 
simulations should provide us with the most reliable information about the temperature 
and/or density dependence of vector meson properties. Currently, these simulations can 
be carried out only for finite temperature systems with zero baryon chemical potential, 
using the quenched approximation. From the hadronic correlation functions at large spatial 
separations the screening masses of hadrons have been obtained in the lattice QCD [T|-^|. In 
a recent paper, Boyd et al || have showed that, up to about 0.92T C , there is no significant 
change of the rho-meson screening mass. They have also found that the quark condensate 
does not change very much until one is extremely close to the critical temperature. This is 
in contrast with the result from chiral perturbation theory which shows a decrease in the 
quark condensate ||. 

Using effective lagrangians that include vector and axial- vector mesons as massive Yang- 
Mills fields of the chiral symmetry, Song || has also found that the rho meson mass at finite 
temperature remains unchanged at the leading order of the temperature. On the other 
hand, Pisarski || has shown using the gauged linear sigma model that the rho meson mass 
increases while the a>i meson mass decreases with temperature, and they become degenerate 
when chiral symmetry is restored. 

To the contrary, based on the restoration of scale invariance of QCD for low momentum 
scales, Brown and Rho and Adami and Brown || have shown that the mass of non- 
strange vector mesons should be reduced in dense matter. This is supported by studies using 
the QCD sum-rule approach PJlOl, the effective hadronic model invoking NN polarization 



TT| |13[, and the quark-meson coupling model [14 



Experimentally, the properties of vector mesons in nuclear medium can be investigated 
by measuring the dileptons produced from heavy ion collisions. Since dileptons are not 
subject to the strong final-state interactions associated with hadronic observables, they are 
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the most promising probe of the properties of hot dense matter formed in the initial stage 
of high energy heavy ion collisions. Dileptons have been proposed as useful observables 
for studying the medium modification of the pion dispersion relation at both finite density 



and temperature [[L^,|16|], the in-medium properties of vector mesons fTT|-pT|l, and the phase 
transition from hadrons to the quark-gluon plasma [P^-|5B|. Indeed, our recent study using 
the relativistic transport model has shown that in heavy ion collisions at SIS/GSI energies 
significant differences exist between the dilepton spectra with and without medium modi- 



fications of the mass and width of vector mesons ||21|| . In particular, we have found, using 
the in-medium rho meson mass from the QCD sum-rule calculation, that the rho peak from 
pion-pion annihilation shifts to around 550 MeV, and its height increases by about a factor 



of four. A calculation similar to that of Ref. [21 1 has been carried out in the relativistic quan- 



tum molecular dynamics (RQMD) ||34|| , and a shift in the rho meson peak in the dilepton 
invariant mass spectrum is also seen in heavy ion collisions at AGS/BNL energies. 

in 



Dileptons have already been measured at Bevalac/LBL by the DLS collaboration p5 
heavy- ion collisions at incident energies around 1 GeV/nucleon. Theoretical studies have 
shown that the observed dileptons with invariant masses above about 450 MeV are mainly 
from pion-pion annihilation |3"E| , |3"T|] . Unfortunately, statistics are not good enough in the 
Bevalac experiments to give definite information on the in-medium vector meson properties. 
However, similar experiments with vastly improved statistics have been planned at SIS/GSI 
by the HADES collaboration |§. 

For heavy ion collisions at the SPS/CERN energies, hot and dense matter is also formed 
in the initial stage of the collisions. We expect that medium effects will lead to a shift 
in vector meson peaks in the dilepton invariant mass spectra from these collisions as well. 
Experiments from both the HELIOS-3 [^] and the CERES j|0] collaboration have shown 
that there is an excess of dileptons over those known and expected sources which cannot be 
explained by uncertainties and errors of the normalization procedures f4"I|| . In particular, 
in the CERES experiment on central S+Au collisions at 200 GeV/nucleon, a significant 
enhancement of dileptons with invariant masses between 250 MeV to 1 GeV over that from 



the proton-nucleus collision has indeed been found. 



In a recent Letter , using the relativistic transport model for the expansion of the fire- 
cylinder formed in the initial stage of S+Au collisions, we have shown that the enhancement 
of low-mass dileptons can be quantitatively explained by the decrease of vector meson masses 
in hot and dense hadronic matter. In this paper we will expand and improve upon Ref. [^2 



We will give the details of our model and the dynamical evolution of the fire-cylinder. The 
7], t]', 0, and a\ mesons that have been omitted in Ref. will be included in our model. 
The medium modification of the rj and rj' mesons are obtained based on their light-quark 
content in free space, neglecting possible medium modifications of their mixing, while a\ 
mesons are introduced based on the Weinberg sum rules, which relate the a\ meson mass 



to that of its chiral partner, the rho meson [43,44]. These additional degrees of freedom are 



useful to reduce the initial temperature and pion chemical potential introduced in Ref. [42 



For a complete comparison with the experimental data for dileptons with invariant mass 
below 1 GeV, we will also evaluate the Dalitz decay of 7r°, t], t]', uj and ai mesons, some 
of which contribute significantly to dileptons with mass below 2m n . We will also study the 
sensitivity of our results with respect to some of the parameters of our model, e.g., the initial 
abundance of mesons and their rapidity distributions. In Section II, we present the extended 
Walecka model in which we couple the light quarks of relevant hadrons to the scalar and 
vector fields in order to provide a consistent way to treat the medium modification of hadrons 
in the relativistic transport model. In Section III, we present the details of the relativistic 
transport model and the dynamical evolution of the fire-cylinder. We will also compare 
our results for particle spectra and rapidity distributions to the experimental data. The 
formalism for evaluating dilepton emission from the expanding fire-cylinder will be given 
in Section IV, together with results, discussions and comparisons with other theoretical 
calculations as well as with the experimental data from both the CERES and the HELIOS-3 
collaboration. Section V is devoted to some general discussions. The paper ends with a brief 
summary in Section VI. 
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II. THE EXTENDED WALECKA MODEL 



To study consistently the effects of dropping vector meson masses on the dilepton spec- 
trum in heavy-ion collisions, we need a model for the in-medium vector meson masses that 
can be incorporated into the relativistic transport model to describe the dynamics of heavy- 
ion collisions. For this purpose we extend the Walecka model from the coupling of 
nucleons to the scalar and vector fields to the coupling of light quarks to the scalar and 
vector fields, using the ideas of the meson-quark coupling model [JTJJ and the constituent 
quark model. Recently, Gelmini and Ritzi jPJ, Brown and Rho |4"7 |, and Furnstahl et al. |48 



have been able to obtain the Walecka mean field model from a chiral Lagrangian. Although 
the latter work appears to be quite different from the former two, the anomalous dimension 



of the loop correction is chosen in |48[ to effectively linearize the Lagrangian as for Gelmini 
and Ritzi j46f and Brown and Rho ^7|. We find that these works give us a firm basis for 
our formulation. 

We consider a system of baryons (we take here the nucleon as an example), pseudoscalar 
mesons (tt and rj mesons), vector mesons (rho and omega mesons), and the axial- vector 
meson (ai) at a temperature T and a baryon density ps- The nucleons couple to both the 
scalar and the vector field, while the nonstrange mesons other than the pion couple only to 
the scalar field, since in the quark model they are made of a quark and an antiquark which 
couple oppositely to the vector field. In the mean-field approximation, the thermodynamical 
potential is then given by j|9 

4 



dk In 



+ 



+ 



+ 



(2tt) 3 
1 

(2^p 
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exp(-(E* N + fi B )/T) + l 
dkln [exp(-£*/T) - 1 
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9 



(2tt) 3 
9 

(2^)3 



exp(-(E* N - n B )/T) + 1 
rfkln \exp(-E n /T) - 1 
exp(-£;/T) - 1 
eM-KjT) - 1 }, 



dk In 
dk In 



(1) 



where (S) and (Vo) are the scalar and vector mean fields, respectively, and the baryon 
chemical potential /x_b is determined by the baryon density 



dtt 4 r r 1 1 -I 

PB ~~d^~Wfl Uxp((££-// B )/T) + 1 ~ exp((^ + ^)/T) + ir U 



In the above expressions, we have -E 1 ^ = y k 2 + m^p and similar expressions for E*, E*, 
E* and E* . Since the medium modification of pions is neglected in this study, we have 



E v = y/W + m*. 

The scalar mean field is determined self-consistently from dQ/d(S) = 0, i.e., 



m ^ s) = §kl d]jr ' 



N 



1 1 



(2n) 3 J E* N Lexp((£ft - (i B )/T) + 1 exp((^ + /x B )/T) + 1 



(2tt) 3 7 £* exp(£>/T) - 1 (2tt) 3 7 E 1 * exp(£;/T) - 1 

29s f„rrC 1 6x/2gg r r< 1 

+ (2vr)3 7 £7* exp(^/T) - 1 + (2vr)3 7 ^ exp(E*jT) - 1' U 

where we have used the constituent quark model relations for nucleon and vector meson 
masses 
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i.e., 

m* N « mjv - <?s(S), m; « m p - (2/3)g s (S), m* « m„ - (2/% 5 (£), (4) 
the quark structure of r\ meson (« 0.58(mt + dd) — 0.57ss) |5(J in free space which leads to 

m*«m„- 0A5g 3 (S), (5) 
and the Weinberg sum rule relation (m ai ~ \^2m p ) between the rho-meson and a\ meson 



masses [E3lE3|, i.e, 



m^«m ai - (2^/3)^(5). (6) 

We have thus neglected in the present model the possible change of rj (also rj') properties at 
high temperature due to the Debye-type screening of instanton effects, which has recently 



been studied in Ref. |pT| , p2| . 

The coupling constants (gv/mv) an d (gs/ m s) are determined by the nuclear matter 
ground state properties. We use here the values corresponding to the original Walecka 
model [fEJ and solve the self-consistent condition to determine baryon and meson masses at 
a given baryon density and temperature. The results for the effective nucleon and rho meson 



masses are shown in Fig. 1. It is seen that for a temperature around 165 MeV and a baryon 
density of about 0.4 fm~ 3 that correspond to the initial conditions of our fire-cylinder, the 
rho meson mass is reduced to about 270 MeV. We would like to emphasize that the effects 
of dropping vector meson masses on the dilepton spectrum do not depend very much on the 
particular model one uses to describe the in-medium vector meson properties, as long as all 
models give similar density and temperature dependence for the vector meson masses. For 
example, Cassing et al. [|53[] have obtained similar results as ours by assuming that the 
vector meson masses in dense matter are given by that from the QCD sum-rule calculations 
@. The choice of the present model is that it offers a consistent way to treat the change of 
vector meson masses during heavy-ion collisions through the change of the scalar potential. 
In this way, the total energy of the system is always conserved. Note that with a stronger 
medium effect in the original Walecka model we have been able to lower the initial baryon 



density substantially from that in Ref. |[42|| , where the nonlinear Walecka model with scalar 



meson self-interactions has been used. The original high baryon density in Ref. gave an 
unrealistically high density of hadrons. 

The above formalism can be generalized to include strange hadrons, such as K, K, and 
hyperons that contain also light quarks. This is achieved by allowing the scalar and vector 
fields to couple to the light quarks in these hadrons. In our model, the masses of these 
particles are thus also reduced at finite temperature and density. Since the phi meson does 
not consist of light quarks, its mass is not modified in the present study. 



III. THE EXPANDING FIRE-CYLINDER 
A. Initial conditions 

Dileptons are produced in all stages of heavy-ion collisions at the SPS/CERN energies. 
Since in this study we are chiefly concerned with dileptons with mass below 1 GeV, the 
initial Drell-Yan processes, that contribute mainly to high-mass dileptons, can be safely 
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neglected. Furthermore, contributions from the quark-gluon-plasma, if it is formed in heavy- 
ion collisions at the SPS/CERN energies, to low-mass dileptons as well as low-momentum 
photons, have been found to be insignificant as compared to those from the hadronic phase 
f5~4] , |55|| . We thus start our transport model from the expansion stage of heavy-ion collisions, 
which might be in the chiral restoration phase as indicated by the decrease of hadron masses. 
Very roughly speaking, this corresponds to the mixed phase in descriptions which do not 
employ medium- dependent masses. The initial hadron abundance and their momentum 
distributions in this expanding fire-cylinder are constrained by requiring that the final proton 
and pion spectra and rapidity distributions reproduce the measured ones. Furthermore, 
results from microscopic simulations based on the string dynamics will also be used to guide 
us in parametrizing the initial conditions [j5"6| -|5"8"|l . We neglect in the present work the isospin 



asymmetry of the colliding nuclei and assume that all charged states of hadrons are equally 
populated. 

1. Baryon abundance 

Let us first discuss the initial baryon density in the fire-cylinder. For a central S+Au 
collision, simple geometrical estimates would give about 90 nucleons from the target nucleus 
if all projectile nucleons are participants. Of course, the number of participant nucleons 
from the target nucleus should be smaller in heavy-ion collisions at the SPS energy of 200 
GeV/nucleon as a result of partial transparency and the not-exactly-zero impact parameter. 
From the experimental data on proton rapidity distributions in the mid- to projectile-rapidity 
region, which have been measured by the NA35 collaboration for S+Au collisions |59| and 



recently by the NA44 collaboration for S+Pb collisions |60| , the total baryon number in 



the fire-cylinder can actually be determined. To describe reasonably the proton rapidity 
distribution, we find that the baryon number in the initial fire-cylinder is about 100 (32 
from the projectile and 68 from the target). The center-of-mass rapidity of the fire-cylinder 
is then 2.65, and this is consistent with the pion rapidity distribution measured by the 
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NA35 collaboration which shows a broad peak around y = 2.6-2.7 We include all 



baryon resonances with masses below 1720 MeV and also the low-lying hyperons, i.e., A, 
E, A(1405) and E(1385). In our model, the hyperons couple to the scalar and vector fields 
with 2/3 of the strength of non-strange baryons, since they contain 2 light quarks instead of 
3 in non-strange baryons. Their initial abundances are determined by assuming strangeness 
saturation as in Refs. ||62|| . Thus knowing the initial baryon density, these abundance can 
be uniquely determined. 

Initially, all baryons are distributed in a fire-cylinder whose cross section is taken to be 
about 40 fm 2 , similar to the geometrical cross section of the projectile nucleus. If we further 
assume that the initial baryon density is about 2.5po (Po = 0.16 fm -3 ), then the initial 
longitudinal length 2zl of the fire-cylinder is found to be 6.2 fm. The initial volume of the 
fire-cylinder is thus about 250 fm -3 . 

2. Meson abundance 

For mesons, we include ir, rj, rj', p, uj, <ft, and ai, as well as K and K*(892). In the 
earlier stage of heavy-ion collisions preceding our initial fire-cylinder, mesons are expected 
to be copiously produced either from string fragmentation as in the RQMD |57|j58[| or from 
the hadronization of a quark-gluon plasma as in ||63|| . They can also be produced from 
processes like nucleon-antinucleon annihilation in a chirally restored phase [Q. The initial 
abundance of (non-strange) mesons (mainly tt, p, to and a\ that feed appreciably to the final 
pion multiplicity) in the fire-cylinder is determined by requiring that the final pion number 
agrees with the measured value and will be specified later. 

For the meson properties in hot and dense hadronic matter, we will use two models. In the 
first model, all mesons are taken to have free masses, and their abundances are determined 
by assuming that they are in thermal and chemical equilibrium, i.e., the 7T, p, u, and a\ 
chemical potentials are related by p p = 2p n and p ai = p w = 3fi n . At a given temperature and 
pion chemical potential, we can determine the meson densities, and from which the meson 



10 



multiplicities are calculated using the initial volume of the fire-cylinder as determined above. 
Normally, one would expect that — so that pions are also in chemical equilibrium with 



baryons, which is indeed the case at the freeze-out as shown in Ref. p2| . However, we find 
that in order to reproduce the observed pion rapidity distribution and spectrum requires 
an initial pion chemical potential of about 130 MeV at an initial temperature of 165 MeV. 
This temperature is determined in our model by fitting the slopes of the proton and pion 
transverse momentum spectra after the full dynamical evolution of the system. Thus, in the 
initial stage of our simulation, mesons are out of chemical equilibrium with baryons, and 
they approach chemical equilibrium as the system expands. Microscopic simulations based 
on the RQMD show that in the central region of heavy-ion collisions the pion density after 
initial string fragmentations is about 0.4 fm -3 , and the pion temperature is about 170 MeV 
f6"5fl . This density is nearly double the equilibrium one and indicates an appreciable pion 
chemical potential for those pions which come from string fragmentation. The initial ir, p, 
u and di meson numbers are then determined to be 119, 56, 43, and 10, respectively. 

Because of finite decay width, the vector meson mass is given the following distribution 

P(M) W , (7) 

1 } (M 2 -m 2 ) 2 + (mr) 2 ' y> 

with a proper normalization. In the above, m denotes the centroid mass of the meson while 
r is its width evaluated at a mass M. 

For rho and phi mesons, their decay widths p — > nn and — > KK are given, respectively, 

by m 

, N g 2 (M 2 - 4m 2 ) 3/2 
W(M) = ^ 12M2 J , (8) 

and 

T^kk(M) - — — . (9) 

where g 2 n7r /A7T pa 2.9 and 9 2 ^ K ^I^ ~ 1-7 are determined from the measured widths at 
m p pa 770 MeV and pa 1020 MeV, respectively. 
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The decay width for K* — > Kir has been derived in Ref. p2| to be 



T K ,^(M) = (10) 

where k is the momentum of the pion in the center-of-mass frame of K*, and g] < „ K7T /A'K w 
0.86 is determined from the measured if* decay width at rriK* ~ 892 MeV. For the decay 
width of ai — > up, we use the result of Ref. i.e., 



T ai ^p = |fe£[2(p, • Pp) 2 + mW + A; 2 )], (11) 



24vrm2 



where k is the pion momentum in the rest frame of a%, and G ai7Tp ~ 14.8 GeV -1 is determined 
from the a\ decay width in free space using its centroid mass. There is, unfortunately, no 
simple expression for the decay width of u — > 7r + 7r~7r° [S3||. We use in this work the 



approximation that this width is directly proportional to the mass of the omega meson, 
which becomes exact in the chiral limit of m n — > 0. 

Since we will consider the r\ and rf Dalitz decay contributions to the dilepton yield, 
their initial multiplicities need to be determined. As for pions, the best way is to fit the 
experimental data after dynamical evolution of the system. There are published data from 
the WA80 collaboration for minimum-biased events and preliminary data for central events 
in S+Au collisions at 200 GeV/nucleon Both are given in terms of 77/71"° ratio 



function of the transverse momentum. We thus determine the initial rj and rj' number mainly 
by fitting the final i]/tt° ratio to the measured one. In the chemical equilibrium scenario 
for the meson chemical composition, both the rj and rf chemical potentials are assumed to 
be about three times that of the pion. The large eta chemical potential is also consistent 



with that of RQMD calculations |65j . This gives us about 37 rj mesons and 5 rf mesons in 
the initial fire- cylinder. The final rj meson number at the freeze-out is about 32 (see Fig. 
6 below). This includes the primary rj mesons that have not been absorbed by nucleons 
or converted into pions through rji] — > tttc, rj mesons from iV*(1535) which decay into rjN 
with a branching ratio of about 0.4 [70], rj mesons from rf which decays into Tnrrj with a 



branching ratio of about 2/3 ffZOfl , and rj mesons produced from nix — > r\r\. We find that the 
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experimental rj/n ratio is reasonably reproduced in this case (see Fig. 10 below). As we 
will show later, in CERES experiments the rj' contribution to dileptons is unimportant but 
that from 77 dominates the dilepton spectrum at invariant mass between 150 < M < 350 
MeV. 

In the second model for the meson properties, we take into account the effects of medium 
modifications according to that predicted by the extended Walecka model as discussed in 
the previous section. At an initial temperature of 165 MeV and an initial baryon density of 
2.5po? we find that the rho and omega meson masses reduce to about 270 MeV, while the 
mass of the a\ meson reduces to about 550 MeV. As a result, their abundance in the initial 
fire-cylinder increases, and we need essentially a zero pion chemical potential in order to 
reproduce the experimental pion rapidity distribution and transverse momentum spectrum. 
We expect a similar trend in RQMD as a part of the energy is now converted to the field 
energy, so less energy is available for pion production, which is, however, compensated by an 
increase of rho, omega, and a x production due to reduced masses. The initial it, p, u, and 
ai meson numbers are then determined to be 44, 92, 33, and 26, respectively. The initial 
low-mass vector and axial-vector mesons thus act as a reservoir for the observed large pion 
yield. As in the free mass case, the mass distributions for these mesons are given by Eq. (|7j) 
with the centroid mass in the equation replaced by the in-medium one. 

In the extended Walecka model, the masses of rj and rj' (using rj' « 0A0(uu + dd) + 0.82ss, 
we have m*, ~ m v i — 0.2(5') [|50|| ) in the initial stage are reduced to about 210 and 810 MeV, 
respectively. Because of the large strange quark content in rj' , the reduction of its mass is 



less significant than that of the rj meson mass as also noted in Ref . [|51[] . When the effective 
masses of rj and r/' are reduced in hot dense matter, either by the screening of instanton 
effects as in Refs. |5T| , |52| or by the attractive scalar field as in our model, their initial 



chemical potentials can be significantly reduced. We are, however, unable to reduce them 
to zero as for the pion chemical potential due to the smaller contribution to rj production 
from higher meson resonances than to pion production. The inclusion of £7a(1) restoration 



5i] , |52|| and ao(980), that decays dominantly into r/7r, will modify this. In addition, since the 
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eta meson mass is reduced in hot dense matter, the process of eta production in pion-pion 
interaction — > r]i] is favored over the reverse one. We find that in order to describe well 
the experimental 77/71"° ratio (see Fig. 10 below), we need initially an eta chemical potential 
of about 130 MeV. The initial 77 and 7/ numbers are then found to be 26 and 3, respectively, 
and the final 77 number at the freeze-out is about 32. 



3. Momentum distributions 

The initial transverse momenta of all hadrons are assumed to be given by a thermal 
distribution, with a proper slope parameter that is loosely referred to as 'temperature' in 
the transverse direction. We find that an initial temperature of about 165 MeV is needed 
to reproduce the observed slopes of the transverse mass spectra for both protons and pions. 



This is somewhat lower than that used in Ref. J42| as a result of the inclusion of additional 
degrees of freedom. Introducing more degrees of freedom should lower it further. The 
longitudinal momentum distribution is determined by imposing a rapidity field as in the 



hydrodynamical model |7T|J72]] . Specifically, we assume that the rapidity of a particle in 
the fire-cylinder frame is correlated to its longitudinal position via the following gaussian 
function, i.e., 

f{y, z) = -jL- exp ( - (y - y L z / z L f /2a 2 ). (12) 
V2iia 



For mesons, RQMD simulations show that their initial rapidity distributions depend 
on the mass, with pions (the lightest) having the broadest rapidity distribution |65|]. For 
simplicity, we will use the same y^ and a in Eq. flT2|) for all mesons. We find that with y^ = 
0.6 we have a good description of the experimental pion rapidity distribution, particularly 
in the mid-rapidity region (see Fig. 8 below), if the variance a is chosen to be 0.8 (also 
for baryons) so that there is substantial dispersion in the y-z correlation. Thus, mesons are 
initially situated mostly around the mid-rapidity but becomes broader at the freeze-out due 
to interactions, leading to good agreements with the experimental data in the mid-rapidity 
region. 
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We note that the initial rapidity distributions in our model is consistent with a momen- 
tum distribution that is in local thermal equilibrium at a temperature of about 165 MeV. 
The introduction of a rapidity field, however, breaks the equilibrium between the longitu- 
dinal and transverse momentum distributions. This is again consistent with that from the 
RQMD 



Baryons, on the other hand, show two components in the initial rapidity distribution [ff3 
The lower component, with jjl = 1.0, is similar to the mesons, while the upper component, 
due to limiting fragmentation, has a larger yi = 2.0. The initial baryon rapidity distribution 
is thus broader than that of mesons. This is necessary in order to fit the observed broad 
proton rapidity distribution, as the proton rapidity distribution is hardly changed by their 
interactions with other particles as a result of their large mass. In Fig. 2, we show the initial 
rapidity and transverse momentum distributions of protons and negative pions for the case 
of free masses. 

B. The relativistic transport model 



The expansion of the fire-cylinder is treated by the relativistic transport model |74|J75|| . In 
this model, hadrons (except for pions whose medium modifications are neglected in present 
work) are propagated in their mean-field potential. For a nucleon, this is 

Jt = W ft=- Vr{E * N+{9V,mv)2pB) ' (13) 

Similar equations are used for baryon resonances with their mean-field potentials taken to 
be the same as the nucleon one. Furthermore, particles also undergo stochastic two-body 
collisions. For baryon-baryon interactions, we include both elastic and inelastic scatter- 
ing for nucleons, A(1232), iV(1440) and iV(1535). Their cross sections are either taken 
from Refs. JSTjlTjJ or obtained using the standard detailed-balance procedure [77]. Lacking 



empirical information, we consider only elastic scattering for interactions involving higher 
baryon resonances and hyperons. This is a reasonable approximation as nucleon, A(1232), 
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iV(1440) and iV(1535) account for about 80% of all baryons. Moreover, in our thermalized 
fire-cylinder, the energy of a colliding baryon pair is not very large, so higher resonance pro- 
duction from baryon-baryon interaction is unimportant. Moreover, their production cross 
sections in baryon-baryon interactions are about an order of magnitude smaller than those 
in meson-baryon interactions, which are included in our model through resonance formation 
and decay. For example, the interaction of a pion with a nucleon proceeds through the 
formation of A(1232), iV(1440), • • •, iV(1720). The formation cross sections are taken to be 
of the relativistic Breit-Wigner form, e.g., the isospin-averaged cross section cr 7r jv- > A(v / s) is 
given by |75| , 



v.n^Ws) - ^ (s _ mi)2 + (mArA)2 > ( 14 ) 

where k is the pion momentum in the center-of-mass frame of the A particle. 

The meson-meson interactions are either formulated by the resonance formation and 
decay when the intermediate meson is explicitly included in our model or treated as a direct 
elastic scattering with a constant cross section estimated from various theoretical models. 

For a pair of pions with a total invariant mass M, a rho meson of this mass is formed 
with an isospin-averaged cross section given by the Breit-Wigner form [21], 



where k is the pion momentum in the center-of-mass frame of rho meson, and the multiplying 
factor (M/rrip) 2 will be clarified later. Similarly, for np — > a± and ttK — > K*, the isospin- 



averaged cross sections 78,79 are, respectively. 



(M\- — (m ai r Q1 ) 2 

and 

(M \_^ (m K *T K *) 2 

k [M 2 — m A K *Y + (m K *T K *y 

Interactions of pions with other mesons are treated through direct elastic scattering. The 

energy-averaged cross section for tci] — > ttt], which is dominated by ao(980) formation and 
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decay, is estimated from the amplitudes given in Ref. |5T| and is about 20 mb, and the same 
is assumed for ttt]' — > tti]'. The energy-averaged cross section for ttuj — > -kuj is estimated 



from the omega collisional broadening width due to the pion-omega interaction in Ref. [80 
and is about 15 mb, and the same is assumed for ttK* — ► nK* and nai — > 7ra\. Finally, we 
have also included the following inelastic processes for pion-pion collisions, tttt — > rji] and 



7T7T — > KK. The cross section for the former is again obtained from Ref. |51| which is about 
10 mb, while the cross section for the latter has been calculated in Ref. [fyj based on the 
one-boson exchange model and is about 1-3 mb depending on the incident energy. 

The energy-averaged cross section for rjr] —>■ rjri is estimated using amplitudes given in 



Ref. [|51[] and is about 12 mb. Interactions of an eta with other mesons are assumed to have 
similar cross sections. The inelastic cross section for 7777 — > titt is obtained from 7T7t — > T]i] 
using detailed balance. 

The collision of a kaon with a pion is treated through the K*(892) formation and decay 
as discussed earlier. The kaon-antikaon collision mainly proceeds through the formation and 
decay of the phi meson, and the isospin-averaged cross section is given by [ F2T 



3tt (m iy) 2 M 2 

The cross section for the direct elastic scattering of a kaon and a rho meson is estimated to 



be about 30 mb from the rho meson collisional broadening width in Ref. ||80fl . This cross 
section is relatively large due to the presence of a broad if* (1270) resonance right above the 
Kp threshold. All other meson-meson collisions are treated as direct elastic processes with 



a constant cross section of 20 mb, which is similar to that used in Ref. [63 



The above discussions apply to both free and medium- dependent masses. In the latter, 
the centroid masses (e.g., m p in Eq. (pT5|)) in these equations are replaced by the in- medium 
ones, while the decay widths (e.g. Eq. (||)) are evaluated with the physical mass M of 
mesons. 
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C. Time evolution of the fire-cylinder 



As an illustration, we show in Fig. 3 the baryon (left panel) and meson (right panel) 
distributions in the x-z plane at several time steps. Initially all particles are located in a 
compressed cylinder whose transverse length is slightly larger than the longitudinal one. 
Then the system expands longitudinally as well as transversely with corresponding flow 
velocities f3 z and (3 t , 



where the summation in the first expression is over all N z particles with longitudinal coor- 
dinates (z — 0.5, z + 0.5) fm, and in the second expression it is over all N t particles with 
transverse coordinates (r t — 0.5, r t + 0.5) fm (r t = y/x 2 + y 2 ). The results are shown in 
Fig. 4 at the freeze-out time (t=20 fm/c). It is seen that the longitudinal flow velocity is 
essentially linear in z. For small r t , the transverse flow velocity is also approximately linear 
in the transverse distance. 

The time evolution of the baryon abundance is shown in Fig. 5. During initial expansion, 
we see an increase of the abundance of higher baryon resonances. This is mainly due to the 
fact that pions are initially out of chemical equilibrium with baryons, and higher baryon 
resonance formation is thus favored over the reverse process. As the system expands, baryon 
resonances eventually decay into pions and nucleons. We do not see much difference between 
the results using free meson masses (the left panel) and in-medium meson masses (the right 
panel) except that the increase of higher baryon resonances is less significant in the latter 
case. This is reasonable as in this case the initial pion density is smaller. A similar plot is 
shown in Fig. 6 for the time evolution of the meson abundance. 

The time evolution of the central baryon density of the fire-cylinder is shown in Fig. 
7, and it is seen that the system expands slower when in-medium meson masses are used. 
In this case, some of the energy goes into the field energy (m 2 s /2){S) 2 . Since this energy 
is independent of baryon number, the associated pressure is negative and thus reduces the 




Pit 



(19) 
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kinetic pressure. This is not a large effect at SPS energies but will be very important at AGS 
energies, where the fraction of energy in the field is appreciable. We plan to work this out 
quantitatively for AGS energies, in the hope that dileptons can be measured at some time in 
the future in this energy domain. The freeze-out temperature of the system, extracted from 
the pion transverse mass spectrum by removing the effect from transverse flow, is about 130 
MeV. The apparent slope parameters of the pion and proton transverse mass spectra are, 
however, larg result of transverse expansion. 

The final proton and pion transverse momentum spectra are compared with the exper- 
imental data in Fig. 8. The results using free meson masses are shown by the dotted 
histograms, while those using in-medium masses are shown by the solid histograms. The 
pion transverse momentum spectra in both cases are seen to be in good agreement with the 
measured data from the CERES |84| and the WA80 |67j collaboration which have been pre- 
sented with arbitrary normalizations. The inverse slope parameters of the pion and proton 
transverse mass spectra have been extracted from our calculation; they are, respectively, 150 
and 220 MeV and are in good agreement with those extracted by the NA44 collaboration 
for central S+Pb collisions at 200 GeV/nucleon |JIJ. 

In Fig. 9, we compare the pion and proton rapidity distributions with the experimental 
data. In our calculation we have not removed the protons from the hyperon decay as in 



the NA44 collaboration [ 60j . We have thus added the hyperons back to the proton rapidity 



distribution given in Ref. by the NA35 collaboration in order to make a meaningful 



comparison. In the mid-rapidity region covered by the CERES collaboration, our pion 



rapidity distribution is in good agreement with the data from the NA35 collaboration pi 
Since we have not attempted to describe the physics around the target rapidity region, which 
does not affect our study of dilepton production in the mid- to projectile-rapidity region, our 
pion rapidity distribution in the small y region should be compared with the open squares 
which are obtained by reflecting the experimental data at large y with respect to y cm ~ 2.65. 
In the mid-rapidity region our pion rapidity distribution also agrees with the preliminary 
data from the NA44 collaboration for central S+Pb collisions |6(| . 
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In Fig. 10, we compare the 77/71"° ratio obtained in our model with the experimental data 
from the WA80 collaboration |68|.|69|]. For the minimum-biased events (solid circles) the 
data have been officially published ||69| , while for central collisions (solid squares) the data 
are very preliminary ||S8|| . Our 77/71"° ratio is larger than the minimum-biased data but 
slightly smaller than the central collision data, and this is especially so in the low transverse 
momentum region. We believe that this corresponds reasonably to the 77/71"° ratio in the 
CERES experiments. With an average charged-particle multiplicity of dN c h/dr) ~ 125 in 
the rapidity range 2.1 < 77 < 2.65, the CERES experiments select a smaller centrality than 
the WA80 experiments which have an average charged-particle multiplicity of more than 



160 in the same rapidity region [^,^]. This gives us the confidence that our results for the 
dilepton production from the 77 and 77' Dalitz decay are quantitatively correct, since their 
contribution to the dilepton yield is essentially determined by their freeze-out multiplicities 
as a result of their much longer lifetime as compared to the lifetime of the fire-cylinder or 
the nucleus- nucleus interaction time (see discussions in next section and Fig. 18 below). 

By fitting and extrapolating the measured 77 and 7r° yields to the full phase space, an 
77/71"° number ratio of 0.147±0.017(stat.)± 0.015(syst.) has been obtained for the minimum- 



biased S+Au collisions by the WA80 collaboration |J§|. This is to be compared with a ratio 
of about 0.08 in the proton-proton interaction at similar center-of-mass energies ||. On the 



other hand, it has been estimated in Ref. 681 that for the 9% most central S+Au collisions 



this ratio could be twice that in minimum-biased collisions, implying an 77/71"° number ratio 
of about 0.29 in the 9% most central collisions. The 77/71"° number ratio in the CERES 
experiments should then be somewhere between these two limits, and our 77/77° number 
ratio of about 0.2 is thus quite reasonable. 

IV. DILEPTON PRODUCTION 

The main contributions to dileptons with mass below 1 GeV are the Dalitz decay of 7r°, 
77 and a), the direct leptonic decay of vector mesons such as p° and u>, and the pion-pion 
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annihilation which proceeds through the p° meson by the vector dominance. In addition, 
we will also evaluate dilepton production from rf and a x decay as well as from decay and 
kaon-antikaon annihilation. 



A. Dalit z decay 

The Dalitz decays of ir°, t], and uo contribute significantly to dileptons with mass below 
2m n . The differential width of the Dalitz decay can be related to its radiative decay width 
87 1 - For example, for ir° — > 7e + e _ , we have |8j 



rfr(7r° -> 7e+e~) _ 4a r(vr -> 27) ( Am 2 A 1/2 



dM 3tt M V M 2 

^IfX'-^V^X (20) 

where M is the mass of the produced dilepton, a is the fine structure constant, r(7r° — > 27) « 



7.65 eV is the neutral pion radiative decay width |70| , and mi is the mass of the lepton. In the 
case of dielectron production, mi = m e « 0.51 MeV can be neglected. The electromagnetic 
form factor is parameterized as 

F n o(M 2 ) = 1 + b w oM 2 , (21) 
with b n o = 5.5 GeV -2 as determined empirically | |87| . For this process the form factor effects 



are small as the dilepton invariant mass is always small (M 
Similarly, for rj / ~fe + e~, we have 

dT(r] —>■ 7e + e") 4a T(r] — > 27) ( Amf ^ 1/2 



x 



dM 3tt M \ M 2 

l + Sf'-SV.lM 1 )! 1 , (22) 



M 2 1 \ m 2 



where T(?7 — > 27) 0.463 KeV is the radiative decay width of an eta meson [7D|. The 



electromagnetic form factor of an eta meson is expressed in the pole approximation as 

F n (M*) = (l-*£)-\ (23) 
v 
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with the cut-off parameter ~ 0.72 GeV as determined experimentally from t] — > / -/p + p^ 
This value is close to that predicted by the vector dominance model and used in Ref. 
In this case the form factor effect is important for dileptons with mass close to m v . 
For <jj — > 7r°e + e~ , the differential decay width is given by 

dT(u -> 7r°e+e-) 2a T(u -> vr°7) , 4mfu/2 / 2m 2 > 



dM 3tt M V m 2 / V m 2 



X 



:i + 



M 2 , 2 . 2m^M , 2 



m 2 , — m 2 m 2 — m 2 



F W7r o(M 2 )| 2 , (24) 



where r(a> — > 7r°7) = 0.717 MeV is the omega meson radiative decay width [70|. The 
electromagnetic form factor is similar to Eq. ( p2]) with A^ replaced by A w . In the vector 
dominance model, the cut-off parameter A w would have a value of about 0.77 GeV, the rho- 
meson mass. However, from the omega muonic decay |87| its value is found to be around 
0.65 GeV, which will be used in the present work. The use of the vector dominance model 
prediction would underestimate the dileptons with mass close to — m n o. 

For 1]' — > / ye + e~ , we have a formula similar to Eq. (P2"|). However, in this case the vector 
meson pole {m p ) occurs in the physical region allowed for the dilepton spectrum (M < m^/), 
and the simple pole approximation is not valid. Instead we use the original vector-dominance 
form factor plj , which has been shown to reproduce the experimental data reasonably well 
for t]' — > / -ffi + fi~ ||87|| . The treatment of the ai Dalitz decay is slightly different. As we 
have in our dynamical model explicitly the processes a\ <-> 7rp and p — > e + e~, the part 
of ai contribution to dileptons which proceeds through the physical p meson as a two-step 
process has already been included. Thus, in evaluating the a\ Dalitz decay (ai — > ne + e~) 
contribution we need not introduce the vector-dominance model form factor. Otherwise, 
there would be double counting. 



B. Vector meson decay into dileptons 



The direct leptonic decay of vector mesons is an important source of dileptons. For the 
dilepton invariant mass region of interest to this work, we consider mainly the leptonic decay 
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of p°, u and <fi mesons. The decay width for p° — > e + e is given by 

r (m _ 9%e 2 M a 2 m p 

r p o_ e+e -(M) _ —- - (25) 

where M 4 in the denominator arises from the virtual photon propagator, and M in the 
numerator comes from the phase space integration. In obtaining the second expression, 
we have used the vector dominance relation g^ = em 2 /g p7T7T |88[| . In the calculation, we 



use, however, r p o^ e + e -(M) = 8.814 x 10 _6 m 4 /M 3 , with the coefficient determined from the 
measured width at M — m p . Similarly, we have T UJ ^ e + e -(M) = 0.767 x 10~ 6 m^/M 3 and 
r^ e + e -(M) = 1.344 x 10~ 6 m^/M 3 for the omega and phi meson dieletron decay widths, 
respectively. 

Combining Eqs. (H), (|15D, and (^), we obtain 

87ra 2 k 

CTn+TT-^pO^e+e-iM) = (T^-^O • T p 0^ e+e -/T p = I F w (M) | 2 , (26) 

with the pion electromagnetic form factor given by 

4 

771 

\ F ^ M )\ 2 = TW2 ^TX7 rvr ( 27 ) 

(M 2 — m z p Y + {m p T P Y 

This is the dilepton production cross section from pion-pion annihilation commonly used 
in the form factor approach []l5|,pTi|36|,^7| . Thus, the dilepton yield from the pion-pion 
annihilation obtained in the dynamical approach with explicit rho meson formation agrees 
with that obtained in the form factor approach if medium effects on the intermediate rho 
meson are neglected [21]]. This explains the multiplying factor (M/m p ) 2 in Eq. ( ]T5|) (and 



a similar factor for kaon-antikaon annihilation, Eq. ([IB])). In Refs. fl21"| , f42"|| we have used 
the rho meson dilepton decay width T p o^ e + e -(M) = 1 / '3a 2 M (47r / g p7Tn ) . In order for the 
relation Eq. ([2H) to hold, we thus need a somewhat complicated multiplying factor for the 
simple Breit-Wigner cross section as shown in Eqs. ( |TTD and fl2?|) of Ref. |2l|. Neglecting the 
medium modification of intermediate vector mesons, the present way of separating a w + n -^ p o 
and r p o^ e + e - and the way used in Refs. PT| , f4"2"| both give the same results as the form factor 
approach for the pion-pion (and similarly kaon-antikaon) contribution. For contributions 
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from the primary rho and omega mesons, different forms of Y p o^^ e + e - as used here and in 
Ref. |4l| do give rise to differences in dilepton yields with masses far away from the rho and 
omega centroid masses. Since low-mass dileptons are dominated by pion-pion annihilation 
and at high masses kaon-antikaon annihilation becomes significant, we find that the results 
of Ref. [fKJ are only slightly affected when the leptonic decay width with an explicit photon 
propagator is used. 

When medium modifications of vector meson masses are included, the rho meson mass 
m p in Eqs. ® and (|2~5|) is replaced by the in- medium mass m* In a static hadronic matter, 



one obtains from the relation cr 7r + 7r -_ >p o_^. e + e -(M) = a^+^-^pO ■ Y p o^ e + e - /Y p again Eqs. ( |26|) 



and ( 2~7D with m p replaced by m*. We note that in obtaining this result the in-medium 
vector dominance relation, g* 7 = em* 2 /g p7T7T , is assumed. With reduced rho meson mass in 
medium, the vector dominance is thus weakened. This effect has also been shown in studies 



based on the hidden gauge theory including loop corrections at finite temperature [89 



Another appreciable medium effect is the collisional broadening of vector meson widths 
in medium. Haglin |5tJ has shown that, at a temperature of 150 MeV, the collisional widths 
are about 50 MeV and 30 MeV for the rho and omega mesons, respectively. However, the 
collisional broadening will not reduce the total dilepton yield. Adding the collisional broad- 
ening width in the denominator of Eq. ( p7|) requires the inclusion of dilepton production 
also from processes responsible for the collisional broadening. The net effect is thus not a 
reduction of the total dilepton yield but only a broadening of the invariant mass spectrum. 
Since the mass resolution in present experiments is comparable to the collisional broadening 
widths, the final results after being corrected by the experimental resolution will be similar 
whether we include the collisional broadening widths or not. We have thus neglected the 
vector meson collisional widths in the present study. We note that the collisional broadening 
widths of vector mesons are reduced when medium-dependent masses are used [pOi . 
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C. Dilepton emission from the expanding fire-cylinder 



In our model, dileptons are emitted not only during the expansion of the fire-cylinder 
but also at the freeze-out. Let us illustrate the way dilepton production is treated in the 
transport model by considering rho meson decay. Denoting the differential multiplicity 
of neutral rho mesons at time t by dN p o(t)/dM, then the differential dilepton production 
probability is given by 

<w e+e - _ ^a*\ t (U)A + y^-i"), (28) 



dM Jo dM p c c y ' dM T P (M) 

where the freeze-out time tf is taken to be 20 fm when there is little interaction among 
the constituents of the fire-cylinder. For the rho meson, whose lifetime is about 1-2 fm/c, 
the contribution from the first term in Eq. (0), i.e, before the freeze-out, is much more 
significant than those from the second term, i.e., after the freeze-out. This remains to be 
so in the case of dropping rho meson mass. Although the rho meson lepton decay width is 
reduced when its mass decreases in medium (see Eq. fl25|)), the increase in the number of low- 
mass rho mesons in medium is more significant, as shown in Fig. 6, leading to an enhanced 
production of low-mass dileptons compared to the free mass case. The in-medium properties 
of rho meson can thus be studied from the shape of its associated dilepton spectrum. 

Similar expressions as Eq. Q2"8| ) are used for dilepton yields from other sources. For 
omega and phi mesons, whose lifetimes are more or less comparable to the lifetime of the fire- 
cylinder, the determination of their in-medium properties from the shape of their associated 
dilepton spectra becomes difficult (see Figs. 13 and 14 below and also Ref. [f21J|). One could, 
however, observe the interesting phenomenon of rho-omega splitting in the dilepton spectra 
if their masses are reduced in hot dense matter as also mentioned in Ref. ||34|| . 

Similarly, the a% Dalitz decay mainly contributes during the evolution of the system 
due to its large width, while the u Dalitz decay has comparable contributions during the 
evolution and at the freeze-out. On the other hand, for the Dalitz decay of 7r°, t], and rj', 
whose lifetimes are much longer than the lifetime of the fire- cylinder, contributions from the 
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first term in Eq. (|28|), i.e., before the freeze-out, are negligibly small (see Fig. 18 below). 
This immediately implies several important observations. First of all, there is almost no 
possibility to observe the medium modification of, e.g., 7] and/or rf meson, from the shapes 
of their associated dilepton spectra, as they decay into dileptons long after the disintegration 
of the hot dense matter, when their properties must have returned to those in free space. 



This has also been realized in Ref. ||52|| . Secondly, the dilepton yield from the Dalitz decay 
of 7T°, T] and r/ is essentially determined by their abundance at the freeze-out. Thus, they 
can be accurately extracted once the abundance of these mesons are reliably determined 
experimentally from, e.g., two-photon reconstruction. Finally, any theoretical calculation 
(or speculation) of their contributions to dilepton production have to be constrained by 
their measured abundance by, e.g., the WA80 collaboration |67|-|69|, as discussed at the end 
of last section. We will come back to this point later. 



D. Results and comparison with the CERES data 

In this subsection, we discuss dilepton production in central S+Au collisions as obtained 
in our calculation. We will show the dilepton invariant mass spectra with and without 
the medium modification of vector meson masses, and with and without the experimental 
acceptance correction. A mass resolution of about 10% of the dilepton invariant mass, in 
accordance with the CERES experiment, has always been included. In Fig. 11, the dilepton 
invariant mass spectra from the Dalitz decay of ir°, r], 00, rj' and a±, as well as the sum of 
these contributions (solid curve), are shown. It is seen that for dileptons with mass below 
100 MeV 7T° Dalitz decay dominates. In the mass region from 150 to 400 MeV, the eta Dalitz 
decay is the most important, while omega Dalitz decay dominates the Dalitz background 
from 500 to 650 MeV. The contributions from rf and a\ Dalitz decay become important 
only above 650 MeV, where the contributions from direct decay of vector mesons are much 
more significant. We do not see much difference between the results with and without the 
medium modification of meson properties. Although more mesons are present because of 
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reduced masses in the case where medium dependence is taken into account, the chemical 
potential enhances their numbers in the free mass case. The dilepton invariant mass spectra 
after applying the experimental acceptance cut are shown in Fig. 12. Note that both the 
rj and uj Dalitz decay contributions in our calculation are somewhat larger than those in 
the 'cocktail' of the CERES collaboration P0| , as we have larger 77/71"° and uj/tt° ratios than 
what were assumed in Ref. f|0| based on proton-proton and proton-nucleus data due to their 



enhanced production in heavy-ion collisions. 

The dilepton invariant mass spectra from the decay of rho (including rho mesons formed 
from pion-pion annihilation, from the decay of a>\ mesons and higher baryon resonances, as 
well as primary rho mesons) and omega mesons are shown in Fig. 13. The dotted curves 
are obtained using free vector meson masses, while the solid curves are obtained using in- 
medium masses. In the low-mass region, the enhancement is about a factor of 3-5 and is 
similar to that found in Ref. |21j for heavy- ion collisions at SIS/GSI energies and in Ref. 



34[ at AGS/BNL energies. This is mainly from the decay of low-mass primary rho mesons 
in the initial hot dense fire-cylinder. We note that as the omega meson has a longer lifetime 
it decays mainly at the freeze-out where its mass has returned to that in free space, leading 
to the peak around m w = 783 MeV in the dilepton invariant mess spectrum. A similar plot 
after the experimental acceptance cut is shown in Fig. 14. 

The sum of dilepton invariant mass spectra from Dalitz decay and direct vector meson 
(including also phi meson) decay from Figs. 11 and 13 (without experimental acceptance 
cut) are shown in Fig. 15. The dotted curve is obtained with free meson masses, while the 
solid curve is obtained with in-medium masses. It is seen that in the mass region from 2m n 
to about 550 MeV, there is an enhancement of about a factor of 3 when the in-medium 
vector meson masses are used. The enhancement at low masses is smaller than shown in 
Figs. 13 and 14 as at these masses the Dalitz decay contribution, which is very similar in 
the two cases, becomes important. 

The comparison with the experimental data is shown in Fig. 16. The results using the 
free meson masses are shown with the dotted curve. It is seen that this model underestimates 
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the data from 250 to 550 MeV by about a factor of 3, while it slightly overestimates the data 
around m PtLU . Most importantly, the shape of the dilepton spectrum is in clear disagreement 
with the experimental data. The results of our calculation using in-medium meson masses 
are shown in the figure by the solid curve. The agreement with the experimental data is 
greatly improved as a result of the factor of 3-4 enhancement in the yield of dileptons with 
masses from 250 MeV to about 550 MeV, and some reduction around m p ^. 



E. Comparison with other models 

In addition to our calculation using the relativistic transport model, several different dy- 
namical models have recently been used to study dilepton production in heavy-ion collisions 



at CERN/SPS energies. In Ref. |55| a hydrodynamical model, with the assumption of quark- 
gluon-plasma formation and a first-order phase transition, has been used for the expansion 
of the fire-cylinder. The results for the dilepton invariant mass spectrum are reproduced 
in Fig. 17 by the dotted curve (the outstanding peak around 1 GeV would disappear if an 
experimental mass resolution of about 100 MeV were properly included in the calculation 
of Ref. ||55||), together with our results using free meson masses (solid curve). In Ref. |]53 



the string dynamics for the initial stage of heavy-ion collisions has been combined with the 



hadronic transport model to study dilepton production. The results of Ref. using free 
meson masses are reproduced in Fig. 17 by the dashed curve. It is interesting to see that 
the dilepton spectra obtained in all three calculations based on different dynamical models 
agree with each other within about 20%, and all underestimate the data at the low- mass 
region by about a factor of 3 and overestimate the data around the rho-meson mass. 

The failure of the three calculations using free vector meson masses and the quantitative 
explanation of the data with in-medium vector meson masses suggest that the observed 
enhancement of low-mass dileptons by the CERES collaboration might be the first direct 
experimental evidence for the decrease of vector-meson masses in hot and dense hadronic 
matter as a result of partial restoration of chiral symmetry. This is indeed a significant 
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conclusion and one may ask whether there are 'conventional' or 'alternative' mechanisms to 
explain the data. 

Conventionally, one might think of the medium modification of the rho meson spectral 
function (or equivalently pion electromagnetic form factor) as a typical nuclear many-body 
effect. It has been shown in Refs. [9~T,52| that the rho-meson spectral function is modified 



in dense nuclear matter as a result of the softened pion dispersion relation due to delta- 



particle— nucleon-hole (AiV -1 ) polarization. The form factor as determined in Ref. [92 
has been used in Ref. J53| to study its effects on the dilepton spectrum. Indeed, one sees 
the enhancement of low-mass dileptons as compared to the case with the free-space rho- 
meson spectral function. A similar result has also been obtained in Ref. |93] in which not 



only the delta-particle— nucleon-hole polarization but also the nucleon-particle— nucleon-hole 
polarization in nuclear matter at finite temperature has been included. However, the effects 



found in Refs. [pq , |93[| , though considerably smaller than the effects from the decreasing rho 
meson mass, are most likely overestimated. Firstly, in heavy-ion collisions at SPS/CERN 
energies, the system is highly excited, so it consists not only of nucleons but also of baryon 
resonances, especially A(1232). The softening of the pion dispersion relation due to the 
delta-particle— nucleon-hole polarization is largely cancelled by the counter-effects from the 
nucleon-particle— delta-hole polarization. In this case the softening of the pion dispersion 
relation depends on an effective density p e ff = Pn — \pa as shown in Ref. f94|. At an 
initial temperature of 165 MeV and initial baryon density of 2.5po, we find that p^ and 
Pa are about 0.8p and 0.9p , respectively, thus the effective density p e ff ~ 0.6p - At this 
density the change of the rho spectral function is already insignificant pl| , |92| . Including 
the counter-effects from higher resonances such as iV*(1440), the effective density will be 
further reduced. Moreover, pions outnumber nucleons by about a factor of 5-10 in heavy- ion 
collisions at CERN/SPS energies, so not all the pions can be converted into (AiV _1 )'s. In 
the calculation of Refs. pl] , p2| , the pion dispersion relation has been obtained always by 



putting only one pion into an infinite nuclear matter. Also, the mismatch in the pion and 
nucleon rapidity distributions in the initial stage when dileptons are produced would reduce 
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the effect of delta-hole polarization. For heavy-ion collisions at SPS/CERN energies, with a 
correct pion dispersion relation as a function of not only nucleon but also baryon resonance 
and pion densities as well as a proper consideration of the nucleon momentum distribution, 
this nuclear many-body effect is expected to be much smaller than that found in Ref. fl53| , |93" 



The pion dispersion relation is also softened at finite temperature due to coupling to 
the pion-rho loop [|95l| . The enhancement of low-mass dileptons from this effect is, however, 
largely offset by the reduction in the pion electromagnetic form factor at finite temperature 

Alternatively, in connection with the interesting possibility of U(1)a symmetry restora- 
tion that leads to medium modifications of the 77 and/or rj' meson masses, several speculations 



have been put forward |51| , [52[| . Looking at the Dalitz decay background, one might suggest 
that the medium modification of the 77 and/or rj' meson masses might enhance low-meson 
dileptons. However, as mentioned previously, the Dalitz decays of the 77 and rj' mesons occur 
dominantly long after the freeze-out of the system. This is quantified in Fig. 18 where 
the dilepton spectrum from the decay of eta mesons up to 20 fm/c is shown by the dotted 
curve, while their total contribution is shown by the solid curve. It is clearly seen that their 
contribution to the dilepton yield before the freeze-out, i.e., in hot and dense medium where 
77 and rj' masses might be modified, is negligible as compared to that after the freeze-out 
when their masses are the same as in free space. Therefore the shape of the dilepton spectra 
from the Dalitz decay of the 77 and rj' mesons remains essentially the same, be there medium 
modifications or not. On the other hand, one might argue that the reduced 77 and 7/ meson 
masses increase their production cross section and thus increase overall the dileptons from 
their Dalitz decay. However, one needs to keep in mind an important constraint of the model 
as provided by the experimental data on the 77/71"° ratio |pT|-|69|j . Since we do reproduce the 
77/71"° ratio reasonably well (see Fig. 9), we believe that our results for the Dalitz decay 
of the 77 and 77' are quantitatively correct. As pointed out earlier, there is already some 
enhancement of Dalitz background with respect to the 'cocktail' of Ref. [|0| . 



In Ref. ||51|1 , partial U(1)a restoration has been studied, and both the 77 and 77' meson 
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masses are found to decrease at high temperature. It has been found in Ref. |5T| that in 
order to describe quantitatively the 77/71° ratio in central S+Au collisions as measured by 



the WA80 collaboration p7|-p9|, one needs a finite 77 chemical potential and a reduced 77 
meson mass. This is in agreement with our findings. We need either a large 77 and rf 
chemical potential, or a moderate chemical potential with reduced 77 and 77' meson masses. 
The enhanced 77 yield in central heavy-ion collisions as compared to that in proton-proton 
and proton-nucleus interaction are, however, not enough to explain the enhanced low-mass 
dilepton yield in the CERES experiments. 

Similarly, U(1)a symmetry restoration in hot hadronic matter has been studied in Ref. 
|52|| , in which the medium modification of rf mesons and their enhanced production are 



emphasized. By increasing the 7/ contribution in the 'cocktail' of Ref. [|(J by 16 times, the 



authors of Ref. |52] found that the Dalitz decays, together with hadron-hadron (especially 



pion-pion) contributions, might be able to explain the enhancement of low-mass dileptons in 
the CERES experiments. It should, however, be emphasized that the 77' contribution in the 
'cocktail' of Ref. j|D| has been obtained with an 77'/ 77° ratio of about 0.04-0.05 [ . Thus a 



16-fold increase of 7/ abundance implies an 77'/ 7r° ratio of about 0.7. Together with a similar 
77/77° ratio, this gives rise to a final 77/77° ratio of greater than 1. This grossly overestimates 
the measured ratio in the most central S+Au collisions by the WA80 collaboration and is 
more so for less central CERES experiments. 

F. Sensitivity of dilepton spectra to initial conditions 

The main uncertainty in our model is the initial conditions for the meson chemical 
composition, the transverse momentum distributions of hadrons, and the rapidity field. 
Microscopic models based on string fragmentations can in principle provide the necessary 
information ||65|| . In the present work, only general features of the results from these mod- 



els have been used. Values of the parameters in the model, particularly the initial meson 
compositions, are determined by fitting the measured proton and pion rapidity and trans- 
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verse momentum distributions. With free vector meson masses, these initial conditions fail 
to explain the observed enhancement of low-mass dileptons. One may wonder if different 
initial meson compositions can be used to increase the dilepton yield without invoking any 
medium effects. Fig. 17 has already demonstrated that using an initial meson composition 
from either string fragmentation or in chemical equilibrium can not lead to an enhanced 
production of low-mass dileptons. Here, we will show the results from a calculation with an 
initial meson chemical composition that is determined by the spin-isospin degeneracy as in 
Ref. [p6fl , i.e., the initial ratio of the n, p, u, and a\ abundance is 1:3:1:3. To reproduce the 
observed pion multiplicity then requires that the numbers are 21, 63, 21, 63 for it, p, u, and 
ai, respectively. The dilepton invariant mass spectrum from this scenario is shown in Fig. 
19 together with the results using the chemical equilibrium scenario. The dilepton spectra 
in both calculations are essentially the same, with the peak around m p ^ slightly lowered 
and a small enhancement around 300 MeV to 600 MeV in the spin-isospin scenario. 

We have also studied the dependence of the dilepton invariant mass spectrum on the 
parameters yi and a, which characterize the initial meson rapidity distributions. By varying 
their values but still requiring that the final pion rapidity distribution agrees reasonably with 
the measured one, we have not found any appreciable change in the dilepton invariant mass 
spectra for both free and in-medium masses. 

G. Comparision with the HELIOS-3 data 

Enhancement of dilepton yields in heavy-ion collisions in both low-mass and 
intermediate-mass regions have also been observed by the HELIOS-3 collaboration which 
have measured dimuon invariant mass spectra in both S+W and p+W collisions at 200 
GeV/nucleon E9]. The enhancement has been demonstrated in terms of the difference in 



dilepton yield per charged particle between S+W and p+W collisions as shown in Ref. ||39 
The data then show that the enhancement is most pronounced around M ~ 0.5 — 0.6 GeV, 
which provides another possible indication, independent of the CERES data, that the vector 
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(rho) meson mass might be reduced in hot and dense matter. 

To see whether our model, that has explained the enhancement of dilepton yield as 
observed by the CERES collaboration, can also explain the enhancement observed by the 
HELIOS-3 collaboration in the low-mass region, we have carried out two calculations for 
dimuon production. One assumes free meson masses and the other uses in-medium masses 
with the same parameters as specified earlier. Since the HELIOS-3 collaboration measures 
dimuons, with m M « 0.106 GeV, we need to include a phase factor of (1 — (2m^/M) 2 ) 1 ^ 2 (l + 
2ml/ M 2 ) in Ec l- (HI) for vector 

meson decays. The results are shown in Fig. 20, and 
compared with the data from Ref. ||39|| . Note that the data are the difference between 
dilepton yields per charged particle in S+W and p+W collisions, i.e., (AN^/N ch ) sw — 
pW> while our theoretical results include only dileptons from the decay of rho and 
phi mesons formed from pion-pion and kaon-antikaon annihilation, which are absent in p+W 
collisions. The y-axis, AN^/Nch (50MeV/c 2 ) -1 , represents the total number of dimuon pairs 
in a mass bin of 50 MeV, including the experimental acceptance cuts and normalized to the 
total charged-particle multiplicity in the pseudorapidity region 3.7 < rj < 5.2. It is seen that 
with free rho meson mass, the shape of the pion-pion contribution is in clear disagreement 
with the data; in the mass region below 650 MeV it underestimates the data by about a 
factor of 2-3, while around m p it slightly overestimates the data. This is very similar to the 
case for the CERES data. Including decreasing rho-meson mass in hot and dense matter, the 
agreement with the data is again greatly improved. Specifically, we have about a factor of 2 
enhancement in the mass region below about 600 MeV as compared with the case using the 
free rho-meson mass. There is also some reduction around m p due to the shift of rho-meson 
peak to a lower mass. Overall, we have reasonably agreement with the experimental data up 
to about 800 MeV. In the mass region around 1 GeV, enhanced dilepton production in S+W 
collisions has also been observed which is most likely due to the kaon-antikaon annihilation 
that proceeds through the phi meson [ f2Tf . 

In both experiments |5S|-^TJ and our theoretical calculation, the enhancement of dilepton 
yield in the low-mass region in Fig. 20 is somewhat smaller than that in Fig. 16. There 

33 



are two possible reasons. Firstly, the CERES collaboration measures dileptons in the mid- 
rapidity region where initially there are many low-mass rho mesons, while the HELIOS-3 
collaboration measures the forward-rapidity region with a smaller number of low-mass rho 
mesons. Secondly, in our model the dileptons measured in the CERES experiments are 
produced early in the expansion when the baryon density is high, and the rho meson mass 
is small. In the HELIOS-3 experiment, the dileptons measured in the forward rapidity are 
produced later in the expansion when interactions lead to a broader rapidity distribution. 
The baryon density at which large rapidity dileptons are produced is thus lower, and the 
reduction of the rho meson mass is smaller. 



V. GENERAL DISCUSSIONS 

There have been discussions of processes which might fill in the low-mass dilepton spec- 
trum. One example is the pionic bremsstrahlung process shown in Fig. 21. This process 
would, indeed, contribute to low-mass dileptons. If the lower left p-meson is produced by a 
7-ray, one can see that this is just the electromagnetic polarizability of the pion, with the 
final 7-ray going into an e + e~ pair. It is known, however, from the Das, Mathur and Okubo 



relation |97j that the p-meson does not enter into the polarizability of the pion, although it 
is most important in the pion electromagnetic form factor. A double-p-seagull enters with 
equal magnitude and opposite sign in order to cancel the process of Fig. 21. This means 



that the p and n couple uniquely through the a\ as intermediate state [pq| , and this has 
been employed in our study. 

For systems with nonzero charge, the ++ sigma meson can couple through 7r + 7r~ to a 
single virtual photon and dileptons |M| . This bears some relation to the known possibility in 
Walecka theory at finite density in nuclei that the sigma can couple to the time component 
of the cu-meson, which can then decay into dileptons. Thus, at finite density (and breaking 
of Lorentz invariance) the sigma cannot decay into a single photon or dilepton pair does not 
hold. Wei don makes the argument that the sigma mass will drop with density and/or 
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temperature, so that dileptons from the sigma decay through a single virtual photon will 
fill up the low-mass region of dileptons. However, as long as the vector mesons are made in 
sufficient number, as we find for the process mi <-> p, it does not matter how many additional 
processes are introduced. Very simply, an equilibrium distribution of vector mesons, with 
appropriate Boltzmann factors, is formed. If additional processes are added to make more 
vector mesons, then the reverse processes will remove it. 

The goodness of our initial conditions, mainly the hadron momentum distributions, was 
found by Cassing et al. ]53| and in the RQMD [|65|| . These authors began from initial 



distributions generated from string dynamics without the assumption of specific functional 
forms for the transverse momentum spectra and rapidity distributions. Yet it can be seen 
from our Fig. 17 (and a similar figure from Ref. ||65|| ) that their results are similar to ours 
and to the hydrodynamical ones of Srivastava et al. ||55|| . Our somewhat larger contribution 
at low masses results from our finite pion and eta chemical potentials. 

Finally, one can question our use of Walecka theory, although many people have used it 
at constituent quark level. It is true that the constituent quark effective masses in Walecka 
theory go to zero only as p — > oo, whereas we believe Nambu— Jona-Lasinio theory to be 
more correct in letting the masses go to zero at finite (and not too high) densities. We have, 
however, started the evolution of our fire-cylinder when the constituent quark has about 
one-third of its on-shell mass, and the difference between Walecka theory and Nambu- Jona- 
Lasinio theory should not be too great here. It can be shown |100|| that the field energy 



in the Walecka theory acts like a bag constant, and that the parameters in the Walecka 
linear sigma-omega model are such that if m*Q is brought to zero, this effective bag constant 
faithfully reproduces that from lattice gauge calculations in QCD. The latter is obtained from 
changes in quark/gluon condensates as the temperature goes through the chiral transition 



101|| . In our expansion of the fire-cylinder it is important to have the correct value for 
the bag constant. The energy going into the bag constant takes away from the energy in 
expansion, providing negative pressure, so the expansion is slower as can be seen from our 
Fig. 6. 
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As noted above, RQMD-type string dynamics calculations have been and are being 
carried out to establish initial conditions. In our dropping mass scenario, at the chiral 
phase transition the vector meson masses are zero. Thus, there is a period between the 
phase transition and the time we begin evolution of the fire-cylinder in which the masses are 
smaller than those in our calculation. The decay of these very low mass vector mesons should 
not be very important, because of the small phase space and because they will chiefly give 
dileptons in the low-mass region, which is dominated by the background of Dalitz decays. 
Work is, however, underway to propagate the masses described in Nambu— Jona-Lasinio 
theory. This theory does bring the constituent quark masses smoothly to zero, and the 
bag constant can be arranged to match the decondensation energy obtained in lattice gauge 
calculations. 



VI. SUMMARY 

In conclusion, using the relativistic transport model for the expansion stage of the fire- 
cylinder formed in the initial stage of heavy-ion collisions at SPS/CERN energies, we have 
studied the enhancement of low-mass dileptons in central S+Au and S+W collisions as 
recently observed by the CERES and the HELIOS-3 collaboration. The initial conditions 
for the fire-cylinder are constrained by microscopic simulations based on string dynamics 
as well as experimentally measured proton and pion rapidity distributions and spectra. 
We have included the Dalitz decay of 7r°, 77, r/, uo and a\ mesons, the direct decay of 
primary p°, u and <fi mesons, and the pion-pion annihilation that proceeds through the p° 
meson, the pion-rho annihilation that proceeds through the a\ meson, and the kaon-antikaon 
annihilation that proceeds through the meson. We have found that the modification of 
vector meson properties, especially the decrease of their mass due to the partial restoration of 
the chiral symmetry, in hot and dense hadronic matter, provides a quantitative explanation 
of the recently observed enhancement of low-mass dileptons by the CERES collaboration in 
central S+Au collisions and by the HELIOS-3 collaboration in central S+W collisions. In 
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comparison with the latter, we have, however, too few phi mesons. It has been shown in Ref. 



102|| that phi meson mass is considerably brought down by the substantial (positive) (ss) 



present in the fire-cylinder. We are unable, as of yet, to include this in a thermodynamically 
consistent way in the transport model. 

To continue this study, we have already started to carry out a calculation [BS| in which 
the initial conditions are completely given by the string dynamics. Although a similar study 
has previously been carried out by Cassing et al. H53"|j56fl , our approach differs from theirs 
in the treatment of vector meson masses in hot dense matter as discussed in Section II. 
Preliminary results show that the results presented in the present paper remain essentially 
unchanged. This more microscopic approach allows us to make predictions regarding low- 
mass dilepton production from collisions between heavier nuclei, where the medium effects 
are expected to be more appreciable. 
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Figure Captions 

Fig. 1: In- medium mass of nucleon (left panel) and rho-meson (right panel) as a function 
of temperature at several baryon densities. 

Fig. 2: Initial proton and pion transverse momentum and rapidity distributions. 

Fig. 3: Baryon (left panel) and meson (right panel) distributions in x-z plane at several 
time steps. 

Fig. 4: Longitudinal (left panel) and transverse (right panel) flow velocities at freeze-out. 
Fig. 5: Time evolution of the nucleon, A(1232) and higher baryon resonance abundance. 
Fig. 6: Time evolution of the meson abundance. 
Fig. 7: Time evolution of the central baryon density. 

Fig. 8: Proton and pion transverse momentum spectra. Dotted and solid histograms are 
obtained from simulations based on the relativistic transport model with free and in- 
medium vector meson masses, respectively. Open and solid circles are experimental 
data from the CERES [84] and the WA80 [^7] collaboration, respectively. 



Fig. 9: Same as Fig. 8 for proton and pion rapidity distributions. The experimental data 
from the NA35 |JT] are given by open circles, while open squares are obtained by 
reflecting the data with respect to y ~ 2.65. Solid circles are the experimental data 
from the NA44 collaboration pU], 



Fig. 10: rj/ir ratio. Open circles and squares are obtained from simulations based on the 
relativistic transport model with free and in-medium vector meson masses, respectively. 
The solid circles and squares are experimental data from the WA80 collaboration for 
the minimum-biased p9 and central |68[ S+Au collisions, respectively. 



Fig. 11: Dilepton invariant mass spectra from Dalitz decay using free masses (left panel) 
and in- medium masses (right panel). 

44 



Fig. 12: Same as Fig. 11 with the experimental acceptance cut. 

Fig. 13: Dilepton invariant mass spectra from the direct decay of rho and omega mesons, 
using free masses (dotted curves) and in-medium masses (solid curves). 

Fig. 14: Same as Fig. 13 with the experimental acceptance cut. 

Fig. 15: Dilepton invariant mass spectra including both Dalitz decay and direct decay 
of vector mesons. The dotted and the solid curve are obtained using the free and 
in-medium masses, respectively. 

Fig. 16: Same as Fig. 15 with the experimental acceptance cut. The experimental data 
from the CERES collaboration |4(| are shown by the solid circles, with the statistical 
errors given by bars. The brackets represent the square root of the quadratic sums of 
the systematic and statistical errors. (We have followed the request of Itzhak Tserruya 
of the CERES collaboration to show the experimental errors in this way instead the 
linear sum used in the original publications [KJ, which were often misinterpreted.) 



Fig. 17: Comparison of different calculations with free meson masses. The solid curve is 
from the present work, while the dotted and dashed curves are from Ref. |55| and Ref. 
|53|| , respectively. 



Fig. 18: Dilepton invariant mass spectra from rj Dalitz decay. The dotted curve gives the 
spectrum up to 20 fm/c, while the solid curve is the total. 

Fig. 19: Dilepton invariant mass spectra obtained using free meson masses. The dotted 
and the solid curve are obtained, respectively, using the spin-isospin scenario and the 
chemical equilibrium scenario for the initial meson composition. 

Fig. 20: Dimuon mass spectra from pion-pion annihilation. The solid and dotted his- 
tograms are obtained with in-medium and free meson masses, respectively. Solid 
circles are the experimental data for the difference in the dimuon yields per charged 
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particle between the S+W and p+W the 
21: Pionic bremsstrahlung. 



collision from the HELIOS-3 collaboration 
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